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ABSTRACT 


A  multicomponent  reliability  system  in  which  each  com¬ 
ponent  is  either  up  (i.e.,  working)  or  down  (i.e,  failed)  in 
accordance  with  an  alternating  renewal  process  is  considered. 
For  arbitrary  structures  the  following  quantities  are  derived.’ 

(i)  The  average  rate  of  system  failure. 

(11)  The  average  uptime  of  the  system.  >  1 

(ill)  The  average  downtime  of  the  system. 

Further  results  are  also  obtained  in  the  special  case  where 
the  system  structure  is  either  series  or  parallel. 
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MULTICOMPONENT  RELIABILITY  SYSTEMS 


Sheldon  M.  Ross 


1.  INTRODUCTION 


Consider  an  n-component  reliability  system  having  the  property  that  at 
any  time  each  of  its  components  is  either  up  (i.e.,  working)  or  down  (i.e., 
being  repaired).  Let 


!1  if  i 
0  othe 


component  is  up  at  time  t 


otherwise 


We  suppose  that  {X^(t),t  ^  0}  ,  i  ■  1,  n  ,  are  Independent  alternating 
renewal  processes.  That  is,  for  each  i  ,  there  exist  random  variables  U*  , 

Dj  •  J  £  1  •  8uch  that 

!1  if  t  <  U* 

0  if  U.1  <  t  <  U?'  +  D* 

1  t  i  i  i 

1  if  uj  +  dJ  <  t  <  u*  +  d*  +  u* 

0  etc. 

The  random  vectors  ,  j  >_  1  ,  are  assumed  to  be  independent  and  identically 


distributed.  Let 


FA(t)  -  p|Uj  <  t| 

Gi(t)  -  p|dJ  <  t} 
00 

U£  -  J  tdF1(t) 

0 

00 

v±  -  J  tdG1  ft ) 
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and  suppose  that  0  <  <  «•  ,  0  <  <  »  ,  F^O)  ■  G^O)  ■  0  . 

We  suppose  that  whether  or  not  the  total  system  Is  working  at  time  t 
depends  only  on  X(t)  ■  (X^(t),  . ..,  XR(t))  .  In  particular,  letting 

(1  if  the  system  is  working  at  time  t 
0  otherwise 

we  assume  that  there  exists  a  nondecreasing  binary  function  <j>  such  that 

X(t)  -  *0^(0,  ....  xn(t)) 

It  follows  from  the  independence  of  components  that  P{X(t)  -  1}  *  E[4>(X(t))] 

is  a  function  of  E[X.(t)],  E[X  (t)]  and  we  write 

l  n 

P(X(t)  -  1}  -  h(E[X.  (t)  ] . E[X  (t)  ] ) 

1  n 

The  function  h  is  called  the  reliability  function  of  the  system.  In  words, 

h(Pp  ...»  pR)  represents  the  probability  that  the  system  will  function  when 
th 

the  i  component  (independently  of  other  components)  functions  with  probability 
Pj  .  Assuming  (as  we  do  throughout)  that  the  distribution  of  U*  +  D*  is  not 
lattice,  i  ■  1,  n  ,  it  follows  from  known  alternating  renewal  process  re¬ 

sults  that 

vi 

E[X, (t) ]  - - -  as  t  ■+  » 

i  +  vi 

and  thus,  from  the  continuity  of  h  (which  is  easily  shown),  we  obtain  that 
(1.1)  Ptt(t)  -  1)  ■»  hi — .  ....  - ^2— 1  =  h(— 2— \ 

\ui+vi  "n  +  v 

a  result  first  noted  by  Esary  and  Proschan  [3j.  The  function  P{X(t)  ■  1}  is 
referred  to  in  the  literature  as  the  availability  at  time  t  ,  and  thus 
Equation  (1.1)  states  that  the  limiting  availability  is 
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2.  REGENERATION  POINTS 


One  unpleasant  feature  of  the  reliability  process  defined  in  Section  1  is 
that  it  does  not  necessarily  have  any  regeneration  points.  That  is,  there 
need  not  be  any  time  points  at  which  the  process  "probabilistically  starts 
over  again."  This  is  unfortunate  as  the  theory  of  regenerative  processes 
and  the  associated  theory  of  renewal  reward  (or  cumulative)  processes  yield 
elegant  results  concerning  not  only  the  existence  of  limits,  but  also  the 
equality  of  different  types  of  limits  such  as  almost  sure  limits  and  limits 
in  expectation.  However,  there  are  certain  special  cases  for  which  the  pro¬ 
cess  is  regenerative.  One  of  these  being  if  for  each  i  *  1,  ...,  n  ,  either 
F^  or  is  a  mixture  of  an  exponential  and  some  other  distribution.  That 
is,  if  for  each  i  ,  either  F^(t)  or  G^(t)  is  of  the  form 


(21) 


»1(l  -  e  ±t:)  +  (1  -  P1)H1(t) 


where  0  <  pt  <  1  ,  and  is  a  distribution  function,  then  the  process  is 

regenerative  with  a  finite  mean  regeneration  cycle*;  the  regeneration  times 
being  those  times  at  which  the  "exponential  parts"  of  each  component  are  in 
effect.  In  fact  since  the  Pj^  in  Equation  (2.1)  can  be  arbitrarily  small 
(though  positive)  it  follows  that  any  set  of  distributions  i  ■  1,  ...,  n) 

can  be  approximated  arbitrarily  closely  by  one  of  the  forms  (2.1).  As  a  result, 
we  shall  assume  that  our  reliability  process  is  a  regenerative  process  with  fi¬ 
nite  mean  regeneration  cycle.  In  fact  since  our  results  will  only  depend  on 
and  v^  ,  i  ■  1,  . . . ,  n  ,  and  as  any  F^  (or  G^)  can  be  arbitrarily 
closely  approximated  by  a  distribution  of  the  form  (2.1)  it  is  intuitively 
obvious  that  our  results  hold  in  the  general  case.  (In  fact,  it  would  seem 
that  appropriate  applications  of  the  ergodic  theorem  yields  the  same  existence 
and  equality  of  limit  results  as  does  the  assumption  of  regeneration  points; 
thus  proving  our  results  in  the  general  case.) 

*The  finite  mean  regeneration  time  follows  from  the  result  that  a  finite  state 
Markov  chain  has  no  null  recurrent  states. 


3.  RATE  OF  SYSTEM  BREAKDOWN 


We  aay  that  component  1  causes  a  breakdown  at  time  t  If 

(i)  X^t")  -  1 

(ii)  X^t)  -  0 

(iii)  X(t")  *  1 

(iv)  X(t)  -  0 

Let  N^(t)  denote  the  number  of  breakdowns  caused  by  i  in  [0,t]  . 
Notation; 

Let  P  -  (P^,  ....  P  )  .  Define 

(1^»P)  *  (Pj^»  •  ••»  ^ *^i+l ’  **** 

(0ltP)  "  (P^»  •••»  •••* 

Proposition  3.1: 

(a) 


(b) 

Proof : 

From  the  theory  of  renewal  reward  processes  it  follows  that,  with 
probability  1, 


With  probability  1, 


N.(t) 

lim  — - —  -  lim 

t-KD  f*®0 


EfN^t  +  a)  -  N1(t)] 


ElNjCt)]  ~ 


>,i  +  v1 


■(V&)  ~ 


“i  +  vi 


N  (t)  E[N  (t)  ] 

lim  ~~  -  lim - -t - 


E[Number  of  breakdowns  caused  by  i  in  a  regeneration 
cycle] /EfTirae  of  a  regeneration  cycle] 
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Define 


I 


i  1  if  the  breakdown  of 
l 0  otherwise 


i  causes  a  system  breakdown 


Now,  again  by  the  theory  of  renewal  reward  processes  it  follows  that,  with 
probability  1, 


(3.1) 


A-  'T  #  •  «  »  A, 

lim - r - ■  lim 

k  k  k 


E[I1  +  ...  +  Ikl 
it 


■  E [Number  of  breakdowns  caused  by  i 
cycle ]/E[ Number  of  breakdowns  of  i 


in  a  regeneration 

in  a  regeneration  cycle] 


However,  from  the  independence  of  components,  it  follows  that 


(3.2) 


E<y  *  "fi'-h)  -  h(0i-z£i)  ae  J  * 


This  is  true  since  the  j  failure  of  i  will  cause  a  system  failure  if  the 
state  vector  x  at  that  time  is  such  that  (Kl^.x)  *  1  .  <KO^,x)  a  0  ,  We 
then  obtain  (3.2)  since  P{<Kl^,X(t))  ■  1,$  (0^,X(t))  ■  0}  -  P{ (l^,X(t))  =  1) 

-  P{iJ>(0i,X(t))  *  1}  .  Thus  from  (3.1)  and  (3.2)  we  see  that,  with  probability, 


T1  +  +  Tk 


h{h-uj)  - 


Now,  since  lim 


Nt(t) 


exists  with  probability  1,  it;  follows  that 


N.  (t)  N.  (time  of  ktk  breakdown  of  i) 

lim  — - —  »  lim  - 

t-*»  c 


Ic*®  time  of  ktk  breakdown  of  i 


lim 

k 


h  +  *  * '  +  h 


time  of  k^  breakdown  of  i 
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h(1i’£j*)  '  h(°fife) 

"i  +  vi 

Part  (b)  follows  from  (a)  since  Blackwell's  Theorem  (of  renewal  theory)  holds 
for  renewal  reward  processes  (see  [  ]).|| 

Letting  N(t)  denote  the  number  of  breakdowns  in  (0,t)  ,  we  obtain, 
from  Proposition  3.1 

Corollary  3.2: 

(a)  With  probability  1, 


11m 


(„t  + 


(b)  E(H(t  +  .)  -  l(t»  -  .  (ut  +  -  h(°i-7±l)] 


Proof; 


(a)  The  nonlattice  assumption  Implies  that  E[N(t)J  ■  e["  £  N.  (t)l  +  o(t) 

Li-1  1  J 

n 

S*»d  so  (a)  follows.  (It  should  be  noted  that  N(t)  i*  £  N.  (t)  since 

1 

more  that)  1  component  can  be  given  credit  for  a  breakdown  at  time  t.) 

(b)  Follows  from  (a)  since  Blackwell's  Theorem  holds  for  renewal  reward 
processes.  1 1 


We  say  that  component  1  causes  an  uptime  at  time  t  if 

(i)  X^t")  -  0 

(ii)  X4(t)  -  1 
(ill)  X(t’)  -  0 
(iv)  X(t)  -  1 


MWMIWrti  in 


Let  N*(t)  denote  the  number  of  uptimes  caused  by  i  in  [0, t ]  .  By  a  proof 
completely  analagous  to  that  of  Proposition  3.1  we  can  prove 


Proposition  3.3: 


(a)  With  probability  1, 


Remark; 

It  is  interesting  to  note  that  the  proportion  of  downtimes  (i.e.,  breakdowns) 
caused  by  component  i  is  equal  to  the  proportion  of  uptimes  that  it  causes. 
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4.  UPTIMES  AND  DOWNTIMES 

Let  U(t)  denote  the  cumulative  amount  of  uptime  in  [0, t ]  .  From  the 
theory  of  renewal  reward  processes  it  follows  that,  with  probability  1, 


-  P{X(t)  -  1}  0  as  t  -*■  ® 

and  thus.,  from  Equation  (1.1), 


(A.l) 


t  ■+•  08 


(with  probability  1). 


The  reliability  system  will  alternate  between  periods  in  which  the  system 

is  up  and  periods  in  which  it  is  down.  Let  us  denote  by  the  length  of 

th  th 

the  i  up  period,  i  ■  1  ,  and  by  ,  the  length  of  the  i  down  period.  It 

should  be  noted  that  the  ,  i  >  1  (and  similarly,  the  D^)  are  neither  inde¬ 
pendent  nor  <*re  they  identically  distributed  random  variables.  However  it  does 
follow  from  results  presented  in  [2]  that 


11m 

rr**° 


lim  E[U.  +  ...  +  U  ]/n 
i  n 


EtuJ 


and 

D..  +  ...  +  D 

lim  — - -  -  lim  Elh,  +  ...  +  D  ]/n  -  E[D  ] 

_ ^  n  1  n  00 


where  (D^)  is  a  random  variable  whose  distribution  is  the  limiting  distri¬ 
bution  of  U  (D  )  . 

n  n 

Moreover  from  Equation  (4.1)  we  note  that 


n 
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5.  SPECIAL  CASES 
5.1  Series  Structure 

The  reliability  system  is  said  to  be  a  series  system  if  in  order  for  it 
to  be  up  it  is  necessary  that  all  of  the  components  be  up.  Hence,  for  a  series 
system 


h(P)  -  n  p. 

i-1  1 

Hence,  from  Corollary  3.2,  we  obtain  that  the  rate  of  breakdowns  of  a  series 
system  is  given  by 


(5.1) 


Unm.  ;  -h-  J  u. 

1  i-1  “i  +  VJ  i-1  1 


while  the  average  length  of  up  and  down  periods  is 


E[U  ] 


n 

l  l/Uj 

i-1  1 


(5.2) 


n 

l  -  n 


BCD J 


1 


..  ,  I  i/u . 

j-i  -j  +  vj  i-i  1 


Remark: 


n!  ’  [|  1/ul] 


If  all  component  uptime  distributions  were  exponential  then  E[U 

for  all  n  .  It  is  Interesting  to  note  that  this  is  also  the  average  of  system 
uptimes  for  arbitrary  distributions.  (This  result  was  also  obtained  in  [ll  where 
the  model  considered  was  series  but  with  the  provision  that  when  the  system  failed 
the  remaining  components  did  not  age  but  were  held  in  a  state  of  suspended 


\ 
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animation  until  the  failed  component  was  repaired.)  Furthermore,  when  the 
uptime  distributions  are  exponential  then  D^.Dj,  ...  are  independent  and 
identically  distributed  and  thus 


*¥ 

E[DJ  -  E [D  ] 
n  00 


n 

1  -  IT 


J-A  UJ  *  VJ 


n 

n 


for  all 


j-1  j  i-1 


n  . 


It  is  interesting  to  note  that  this  result  (in  the  exponential  component  uptimes 
case)  seems  to  be  difficult  to  prove  directly. | f 


Additional  insight  into  the  length  of  uptime  periods  is  obtained  by  the 
following  heuristic  reasoning.  Given  that  (the  fixing  of)  component  1  has 
just  caused  an  uptime  period  it  is  intuitive  that  the  remaining  uptime  of  com¬ 
ponent  j  ,  j  i  i  ,  should  be  given  by  the  equilibrium  distribution  F  (x)  , 

j  »e 

where 


(x))dx 


F.  (x)  - 

J  i*  U, 


V 


Hence,  as  the  proportion  of  uptimes  that  are  caused  by  component  i  equals 


lim 

t-H» 


wj(t) 


n  *  n 

I  N(t)  l  1/u, 

3-1  1  1  3 


,  it  follows  upon  conditioning  that 


(5.3) 


F{U  >  x)  -  (  l  1/wA  1  I  ^-(1  -  F.  (x))  n  (1  -  F  (x» 
\j-i  J/  i-i  yi  1  ji*i  j,e 


represents  the  limiting  distribution  of  .  Another  interpretation  of  (5.3) 

is  that  it  represents  the  proportion  of  uptime  periods  that  are  of  length  greater 

00 

than  x  .  If  Equation  (5.3)  is  Indeed  valid  then,  as  E[Ul  -  f  P{U  >  x}dx  , 

w  r  00  w 

0 
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Jn  1 

T  —  (1  -  F. (x))  II  (1  -  F.  (x))dx  to  equal  1 
0  i-1  Wi  1  j*  3,6 


We  now  verify  this. 


Proposition  5.1; 


J 


n  1  -  F  (x) 

l  - ^ -  n  (1  -  F  (x) )dx  -  1  . 

i-1  Wi  J^i  3*e 


Proof: 


The  above  appears  difficult  to  verify  analytically,  but  there  is  a  simple 

1  -  F^x) 

probabilistic  proof.  Noting  that  - - - dx  -  dF^  g(x)  we  see  that  the  above 

integral  is  equal  to 


°5» 

l  [ 

i-1 1 


n  (l 
1+i 


F  (x))dF.  (x)  . 

J.e  i.e 


Letting  X. ,  . . . ,  X  be  independent  random  variables,  with  X.  having  distri- 
I  n  l 

but  ion  g(x)  we  note  by  a  simple  conditioning  argument  that 


P{XA  -  smallest  of 


(xr  ....  xn)>  -  J 


R  (1  -  F.  (x))dF.  <x> 

Jl*i  j,e  i,e 


and  thus,  summing  over  all  i  fields  the  desired  result. 


Remark: 


It  also  f 0II0W8  from  Equation  (5.3)  that  if  each  F^  has  a  DFR  distribu¬ 
tion;  that  is,  if  t  t  for  all  s  ,  then  the  distribution  of 


is  also  DFR.  This  follows  from  the  3  facts 


l 


\ 


X 


(!)  Mixtures  of  DFR  distributions  are  themselves  DFR. 
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(11)  If  F  Is  DFR  then  F  Is  slso  DFR  (this  follows  from  (1)  using 

e 

results  about  the  equilibrium  renewal  process,  It  can  also  be 
shown  directly). 

(Ill)  A  series  system  made  up  of  Independent  DFR  component  lifetimes,  has 
a  DFR  lifetime  distribution. 


In  fact  the  same  argument  shows  that  UR  is  DFR  for  all  n  . 

Another  quantity  of  Interest  Is  the  remaining  uptime  of  a  system  that  is 
up  at  time  t  (t  large).  Again,  heuristic  reasoning  yields  that 


(5.4) 


P{addltlonal  uptime  >  x]  ■  II  (1  -  F  (x)) 

j-1  J'e 


and  thus 


E[additlonal  uptime] 


f  n  (l  -  f 

Jn  J-1  J’ 


(x))dx 


1  -  F.  (x) 


,  fili 

“i. 


n  (1  -  F.  (x))dx 


e  J*1 


00 


where  u 


f  f.Vw 

l,e  'J0  xdFl,e(l°  'l  ■  Hence> 


Efadditional  uptime]  -  y.  P{X.  -  smallest  of  (X.  ,X,,  ,  ....  X  )} 

i|6|6  1|6|G  4^6  Ilg  6 


where  X.  ,  X.  ,  1  >  2  ,  are  independent  and  X,  ~  F,  and 
l,e,e  i,e*  —  r  l»c,e  l,e,e 

Xle~F1  e,  1  >  2  .  In  the  special  case  where  all  the  component  uptime 
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distributions  are  Identical  and  the  distribution  is  IFR  j^that  is, 

—  4.  t  for  an  fl|  it  follows  that  the  expected  additional  system 

uptime  of  a  series  system  that  is  up  at  time  infinity  satisfies  the  following 
inequality 


E[addltional  uptime] 


This  inequality  is  reversed  if  F  is  DFR.  The  above  inequality  follows  from 
the  fact  (proven  below)  that  if  F  is  IFR  then  so  is  Fg  . 

It  again  follows  that  if  each  F^  is  DFR  then  the  distribution  of  addi¬ 
tional  life  as  given  by  Equation  (5.4)  is  also  DFR.  Moreover  if  each  F^^  is 
IFR  then  it  follows  that  the  additional  life  distribution  is  also  IFR.  This 
follows  from  the  fact  that  a  series  system  of  independent  IFR  components  itself 
has  an  IFR  life  distribution  upon  application  of  the  following  proposition. 


Proposition  5.2; 


If  F  is  an  IFR  distribution  with  finite  mean,  then  Fg  is  also  IFR  where 


■1 


F  (t)  -  I  (1  -  F(y))dy/  I  (1  -  F(y))dy 


/ 


Proof : 


The  failure  rate  function  of  Fg  is  given  by 


1  -  F  (t) 
x  (t)  -  I- - Z— 

' dt  Fe(t) 


1-1 


•[I 


1.  -  Hy)  dy 

1  -  F(t)  ay 


-1 
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-  1/E [additional  life  of  a  t  year  old  item 
having  life  distribution  F] 

and  the  result  follows  since  F  IFR  implies  that  the  denominator  of  the  above 
term  is  decreasing  in  t  .  |  | 

It  should  be  noted  that  it  has  been  shown  in  [A]  that  if  all  component 
uptime  and  downtime  distributions  are  exponential  then  both  the  limiting  dis¬ 
tribution  of  system  uptime  and  of  the  additional  uptime  of  a  system  that 

is  up  are  both  mixtures  of  exponential  distributions  and  are  thus  DFR.  These 
results  hold  for  Markovian  (i.e.,  exponential)  systems  having  an  arbitrary 
structure  (that  is,  they  need  not  be  series). 

5.2  Parallel  Structure 

The  reliability  system  is  said  to  be  a  parallel  system  if  it  is  up  whenever 
at  least  one  of  its  components  is  up.  That  is. 


n 

h(P)  -  1  -  n  (1  -  P  )  . 

1 

In  other  words  it  is  down  if  and  only  if  all  of  its  components  are  drum.  Thus 
we  see  that  by  regarding  down  as  up  the  parallel  system  is  transformed  into  a 
series  system  and  so  all  of  the  results  of  Section  5.1  are  immediately  trans¬ 
latable  into  results  about  parallel  systems.  For  example,  from  Equations  (5.1) 
and  (5.2)  we  see  that,  in  the  parallel  case, 
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